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INTRODUCTION 

Approximately  75%  of  breast  cancers  express  the  hormone  estrogen  receptor  a  (ER).  As  a  critical  determinant 
in  estrogen  response  and  oncogenic  driver  for  ER-positive  breast  cancer,  ER  promotes  cancer  cell  proliferation, 
survival,  and  metastasis.  Successful  endocrine  therapy  targets  the  ER  pathway  by  inhibiting  estrogen  synthesis 
with  aromatase  inhibitors,  blocking  ER  with  selective  ER  modulators  (such  as  tamoxifen),  or  eliminating  ER 
with  selective  ER  degraders  (such  as  fulvestrant).  However,  in  ER-positive  (+)  patients  with  metastatic  disease, 
more  than  50%  of  patients  fail  to  respond  to  first-line  therapy  due  to  de  novo  resistance,  and  all  patients  who  do 
respond  eventually  relapse  and  die  due  to  acquired  resistance.  The  mechanism  of  endocrine  resistance  has  been 
poorly  understood.  Our  subject  in  this  study  is  endocrine-resistant  (Endo-R)  ER+  breast  cancer  disease. 
Through  integrated  approach  built  on  multi-omics  platforms  including  next-generation  sequencing  of  a  large 
panel  of  our  established  endocrine-resistant  cell  models,  we  previously  identified  a  novel  ER/FOXA1/IL-8 
signaling  axis  in  our  Endo-R  cell  models.  By  completing  this  DoD  funded  study,  we  hope  to  better  understand 
the  mechanism  of  endocrine  resistance  driven  by  the  alterations  of  the  ER/FOXA1/IL-8  axis,  and  further 
develop  novel  therapeutic  approach  to  target  this  axis  to  overcome  endocrine  resistance  and  improve  patient 
outcome.  The  scope  of  this  study  covers  molecular,  cellular,  and  tumor  biology,  and  integrated  bioinformatics 
analysis  of  multi-omics  data  generated  from  both  in  vitro  cell  line  and  in  vivo  xenograft  mouse  models. 

KEYWORDS 

FOXA1,  gene  amplification,  estrogen  receptor,  interlukin-8,  cytokine,  transcriptional  reprogramming,  endocrine 
resistance,  xenograft  tumor,  growth  factor  receptor,  RNA-Seq,  ChIP-Seq,  multi-omics,  integrative  analysis 

ACCOMPLISHMENTS 

1.  Specific  Aim  2:  Determine  the  role  and  therapeutic  potential  of  FOXA1/IL-8  in  endocrine  resistance,  and 
the  mechanisms  by  which  the  ER/FOXA1/IL-8  axis  mediates  endocrine  resistance  and  tumor  progression. 

Major  Task  4:  Determine  if  overexpression  of  FOXA1  confers  resistance  and  if  downregulation  of 
FOXA1/IL-8  enhances  sensitivity  to  endocrine  therapy  in  our  4-5  prioritized  Endo-R  cell  models.  (Months 
13-23) 

Status:  We  completed  100%  of  the  SOW  of  the  2nd  year.  The  work  we  have  done  in  the  2nd  year  was 
summarized  and  published  in  a  recent  PNAS  paper  (Fu,  et  al.,  Oct,  2016).  We  report  here  our  project 
progression  mainly  addressing  the  specific  SOW.  Most  of  these  data  were  published  in  this  PNAS  paper 
attached  in  the  end  of  this  report.  We  reached  the  milestone  —  We  determined  the  cause-and-effect  of  altered 
FOXA1/IL-8  expression  on  endocrine  sensitivity  across  our  multiple  Endo-R  cell  models. 

There  are  no  significant  changes  in  approach  or  methods  from  the  approved  SOW. 

2.1  Subtask  1.  Examine  whether  FOXA1  upregulation  decreases  endocrine  sensitivity  in  the  stable  P  cells  with 
inducible  FOXA1  over  expression  ( constructed  in  Major  Task  2)  by  using  cell  growth,  proliferation, 
apoptosis,  migration,  and  invasion  assays  on  the  FFTP  Celigo  cell  cytometer  platform.  Months  13-19 

We  have  established  Dox-inducible  FOXA1  overexpression  (OE)  in  MCF7L,  ZR75-1,  and  600MPE  ER+ 
BC  cell  line  models.  The  induction  of  FOXA1  upon  Dox  was  confirmed  by  Western  blot  in  all  three  cell 
models.  The  dose-dependent  FOXA1-OE  upon  Dox  was  confirmed  by  Western  blot  in  MCF7L-P/FOXA1 
cells.  After  pre-starvation  in  5%  charcoal-stripped  FBS  and  phenol-red  free  medium  for  5  days,  the  cells 
with  -/+Dox  were  treated  by  estrogen  (E2,  1  nM)  as  positive  control,  or  anti-estrogens  including  estrogen- 
deprivation  (ED),  tamoxifen  (Tam,  100  nM),  and  fulvestrant  (Ful,  100  nM)  for  an  average  7  days.  Cell 
growth  was  measured  by  either  the  HTP  Celigo  cell  cytometer  or  the  colorimetric  Methylene  Blue  staining. 
Cell  growth  at  endpoint  was  subtracted  by  day  0  and  normalized  to  E2  control.  As  shown  in  Fig.  1A  and  B, 
FOXA1-OE  significantly  reduced  sensitivity  to  all  anti-estrogens  in  both  MCF7L  and  ZR75-1  cell  models. 
We  also  performed  the  cell  invasion  assay  in  a  24-transwell  platform  using  the  inserts  coated  with  Metrigel. 
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FOXA1-OE  significantly  increases  cell  invasiveness  in  both  MCF7L  and  600MPE  cell  models  (Fig.  1C). 
We  also  observed  the  same  results  in  migration  assay. 
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Fig.  1.  FOXA1  overexpression  reduces  endocrine  sensitivity  and  increases  cell  invasiveness  in  ER+  BC  cell  line 
models.  A  and  B,  Cell  growth  assay  in  MCF7L-P/FOXA1  and  ZR75-1-P/FOXA1  cell  models  under  endocrine  treatment. 
C,  Cell  invasion  assay  in  MCF7L,  600MPE  and  ZR75-1  cells  with  FOXA1-OE. 


2.2  Subtask  2.  Determine  whether  knocking  down  ER/FOXA1/IL-8  by  the  pINDUCER  system,  or 
pharmacological  IL-8  suppression,  can  overcome  resistance  in  Endo-R  cell  models.  Months  15-21 

We  performed  ER  and  FOXA1  KD  using  siRNAs  that 
were  validated  by  qRT-PCR  and  Western  blots.  As 
shown  in  Fig.  2A-D,  ER-KD  significantly  suppresses  cell 
growth  in  both  P  and  Endo-R  cells  of  MCF7L  and 
MCF7RN  models.  The  less  suppression  in  Endo-R  cells 
of  ZR75-1  and  600MPE  models  compared  to  the  P  cells 
by  ER-KD  may  be  explained  by  the  loss  of  ER  expression 
in  these  two  Endo-R  cell  models.  FOXA1-KD  suppresses 
the  cell  growth  in  all  tested  Endo-R  cell  models,  though 
there  is  less  effect  in  ZR75-1  Endo-R  cells.  These  data 
suggest  that  FOXA1  is  critical  for  both  P  and  Endo-R  cell 
growth  and  in  cells  with  or  without  ER  expression,  thus 
potentially  involving  both  ER-dependent  and  independent 
pathways  for  cell  growth.  We  have  constructed  Dox- 
inducible  FOXA1  shRNA  in  both  MCF7L-P  and  MCF7L- 
TamR  cells,  based  on  one  of  the  siRNA  sequences  that 
was  validated  before.  We  observed  similarly  potent  cell 
growth  inhibition  upon  shRNA  induction  by  Dox  in  both 
P  and  TamR  cells. 

Since  IL-8  is  highly  upregulated  in  MCF7L-TamR  vs. 
P  cells,  we  performed  IL-8  KD  using  two  different 
siRNAs  in  these  cells,  and  observed  significant  cell 
growth  inhibition  in  TamR  vs.  P  cells  (Fig.  2E). 
Importantly,  this  cell  growth  inhibition  in  TamR  cells  can 
be  rescued  by  simultaneous  induction  of  IL-8  expression 
in  TamR  cells  (Fig.  2F),  supporting  a  non-off-target  effect 
of  IL-8  KD  in  these  cells.  We  further  constructed  the  Dox- 
inducible  EL-8  shRNA  KD  in  both  MCF7L-P  and  TamR 
cells,  and  observed  the  same  effect  in  cell  growth  inhibition  especially  in  TamR  cells  upon  Dox  induction. 
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Fig.  2.  ER/FOXA1/IL-8  KD  inhibits  cell  growth  in 
multiple  Endo-R  cells.  A-D ,  Cell  growth  in  6  days 
after  ER/FOXA1  or  non-specific  (N.S.)  KD  in 
MCF7L,  MCF7RN,  ZR75-1,  and  600MPE  Endo-R 
cell  models.  All  normalized  to  N.S.  KD.  E ,  Cell 
growth  in  MCF7L-P  and  TamR  cells  with  IL-8  KD 
by  two  different  sequences.  F ,  IL-8  KD  rescue  assay 
in  MCF7L-TamR  cells  with  inducible  IL-8 
expression.  The  #2  siRNA  targets  the  3’-UTR  of  IL8. 
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These  data  suggest  that  targeting  the  ER/FOXA1/IL-8  axis  may  be  an  effective  approach  to  overcome 
endocrine  resistance  in  ER+  BC. 

Since  cell  growth  inhibition  by  IL-8  KD  in  TamR  cells  can  be  rescued  by 
overexpression  of  IL-8  cDNA,  we  asked  whether  it  also  can  be  rescued  by 
exogenous  IL-8  protein  supplemented  in  the  culture  medium.  However,  no 
rescuing  effect  was  observed  in  our  MCF7L-TamR  cells  treated  by 
simultaneous  addition  of  IL-8  in  culture  medium  upon  IL-8  KD  (Fig.  3A). 

In  line  with  this,  no  cell  growth  inhibition  was  seen  in  both  P  and  TamR 
cells  upon  addition  of  the  neutralizing  IL-8  antibody  (Fig.  3B).  We  also 
tested  two  of  the  neutralizing  antibodies  of  the  IL-8  receptor  CXCR1  and 
CXCR2.  They  also  cannot  inhibit  both  MCF7L-P  and  TamR  cell  growth. 

These  data  suggest  that  the  effect  of  IL-8  on  the  cell  growth  especially  in  the 
TamR  cells  may  be  involved  in  the  intracellular  IL-8  signaling.  In  line  with 
this  hypothesis,  the  IL-8  protein  was  shown  to  be  localized  mainly  in  BC  cell 
cytoplasm,  especially  in  the  para-nuclear  region  (see  attached  PNAS  paper) 
in  our  in  vitro  cell  lines,  in  vivo  xenograft  tumors,  and  clinical  human 
specimens.  It  has  also  been  reported  that  IL-8  together  with  CXCR1/2 
proteins  show  nonapical  and  cytoplasmic  expression  in  prostate  cancer  cells 
and  induce  autocrine  signaling  associated  with  cancer  progression  and 
angiogenesis  ( Clin  Cancer  Res.  PMID:  15930347).  We  also  cannot  exclude 
the  possibility  of  paracrine  signaling  via  the  secreted  IL-8  between  the 
cancer  and  stromal  cells  in  vivo,  which  is  shown  in  the  PTEN-deficient 
prostate  cancer  cells  that  secret  IL-8  to  induce  chemokines  CCL2  and 
CXCL12  from  stromal  cells  for  cell  invasion  in  a  co-culture  cell  system  ( Oncotarget ,  PMID:  24970800). 
However,  to  mechanistically  test  our  hypothesis  of  EL-8  intracellular  signaling  in  Endo-R  BC  is  out  of  the 
scope  of  this  project.  We  plan  to  extend  this  study  by  pursuing  other  research  funding. 

2.3  Subtask  3.  Determine  whether  IL-8  expression  can  rescue  the  effect  of  FOXA1  on  endocrine  sensitivity  by 
using  siRNA  or  inducible  overexpression  in  our  3  prioritized  Endo-R  cell  models.  Months  20-23 

IL-8  expression  is  highly 
induced  by  FOXA1  in 
MCF7L-P  cells,  and  also 
highly  increased  in 
MCF7L-TamR  vs.  P 
cells.  FOXA1-OE 

reduces  endocrine 

sensitivity  in  MCF7L-P 
cells,  partially  due  to  the 
highly  induced  IL-8 
expression.  To  address 
the  role  of  IL-8  in 
mediating  FOXA1  effect, 
we  constructed  double 
lentiviral  Dox-inducible 
systems  in  MCF7L-P 
cells  with  simultaneous 
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Fig.  4.  Reduced  endocrine  sensitivity  and  proteomic  alteration  upon  FOXA1-OE  can 
be  partly  reversed  by  IL-8  KD  in  MCF7L-P  cells.  A  and  B ,  Cell  growth  in  7  days  in  four 
MCF7L  derivative  lines  -/+  Dox  induction,  treated  with  E2  (as  control)  or  anti-estrogen  (ED, 
Tam,  or  Ful).  Luc  and  YFP  are  negative  control  for  KD  and  OE,  respectively.  C,  Heatmap 
of  partial  RPPA  data  of  the  MCF7L/FOXA1  +Dox  cells  with  siRNA  knockdown  of  N.S.  or 
IL-8.  Proteins  were  chosen  based  on  the  previous  RPPA  data  showing  the  high  expression 
levels  in  MCF7L/FOXA1  +Dox  vs.  -Dox  cells  (log2  ratio  >  0,  right  panel). 
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Fig.  3.  The  effect  of  IL-8  KD 
on  cell  growth  inhibition 
may  be  involved  in  an 
intracellular  IL-8  signaling. 

A,  Exogenous  IL-8  protein 
cannot  rescue  cell  growth 
inhibition  upon  IL-8  KD  in 
TamR  cells.  B ,  Neutralizing 
IL-8  antibody  cannot  inhibit 
cell  growth. 
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FOXA1-OE  and  IL-8-KD,  or  either  one  of  them  with  relative  negative  controls  for  OE  (YFP)  and  KD  (sh- 
luciferase,  shLuc).  Upon  -/+Dox  induction  and  pre-starvation  in  charcoal-stripped  FBS  and  phenol-red  free 
medium,  cells  were  subjected  to  E2  treatment  or  anti-estrogens  (ED,  Tam,  or  Ful)  for  cell  growth  assay.  As 
expected,  without  Dox  induction,  endocrine  therapy  suppressed  cell  growth  in  all  groups  without  difference 
(Fig.  4A).  Upon  Dox  induction,  FOXA1-OE  reduced  endocrine  sensitivity  to  all  anti-estrogens,  and  IL-8 
KD  sensitized  cells  to  anti-estrogens  especially  ED  (Fig.  4B).  Simultaneous  IL-8  KD  in  the  context  of 
FOXA1-OE  significantly  re-sensitized  cells  to  all  anti-estrogens,  supporting  the  role  of  IL-8  in  mediating 
endocrine  resistance  upon  FOXA1-OE.  To  further  understand  the  molecular  changes  underlying  the 
observance  of  endocrine  response,  we  analyzed  the  proteomic  profiles  in  MCF7L-P  cells  with  FOXA1-OE 
compared  to  the  cells  with  simultaneous  IL-8  KD  using  RPPA  assay.  We  specifically  focused  on  the 
proteins  that  were  highly  increased  upon  FOXA1-OE  in  P  cells.  Upon  simultaneous  IL-8  KD,  reduced 
expression  was  found  in  over  50%  of  these  proteins,  which  include  many  phosphorylated  proteins 
representing  activated  signaling  in  MAPK,  Akt,  Stat3,  and  mTOR  pathways  (Fig.  4C).  These  data  suggest 
that  FOXA1-OE  induces  growth  factor  receptor  downstream  signaling  and  endocrine  resistance  at  least 
partially  through  inducing  IL-8. 

We  also  examined  the  IL-8 
effect  on  both  cell  growth  and 
cell  invasion  in  600MPE  Endo-R 
cell  model.  There  is  no  obvious 
cell  growth  inhibition  by  IL-8  KD 
in  both  600MPE-P  and  TamR 
cells  (Fig.  5A),  possible  due  to 
the  endogenous  activated  MAPK 
pathway  associated  with  a  KRAS 
gene  mutation  in  this  line.  In 
contrast,  IL-8  KD  significantly 
suppressed  cell  invasion  of  both  P 
and  TamR  cells  (Fig.  5B).  Importantly,  the  increased  cell  invasiveness  upon  FOXA1-OE  was  mitigated  by 
simultaneous  IL-8  KD  in  600MPE-P  cells  (Fig.  5C),  suggesting  that  the  IL-8  expression  mainly  affects  the 
cell  invasiveness  induced  by  FOXA1  in  this  cell  model.  These  data  also  suggest  that  pleiotropic  effects  of 
IL-8  regulated  by  high  FOXA1  may  contribute  to  different  cell  behaviors  of  endocrine  response  and/or  cell 
invasiveness  due  to  the  heterogeneity  such  as  different  genetic  background  of  ER+  breast  cancer. 

3.  Opportunities  for  training  and  professional  development 

Nothing  to  Report. 

4.  Results  disseminated  to  communities  of  interest 

Nothing  to  Report. 

5.  Plan  to  do  during  the  next  reporting  period  to  accomplish  the  goals 

In  the  next  funding  year,  we  will  perform  in  vivo  experiment  to  determine  the  effect  of  targeting  FOXA1 
and/or  IL-8  on  endocrine  response  in  our  ER+  breast  cancer  xenograft  mouse  model.  We  will  depict  the 
signaling  changes  upon  FOXA1/IL-8  alterations  in  Endo-R  cell  models  by  using  RPPA  and  bioinformatics 
analysis.  Finally,  we  will  determine  the  clinical  relevance  of  the  FOXA1/ER/IL-8  axis  by  measuring  their 
expression  in  collected  human  specimens  from  our  collaborators.  All  these  experiments  will  be  performed 
according  to  the  approved  SOW. 

IMPACT 

Nothing  to  Report. 
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Fig.  5.  IL-8  KD  suppresses  cell  invasion  in  600MPE  cells  with  high  FOXA1. 

A,  IL-8  KD  has  no  effect  on  cell  growth  in  both  P  and  TamR  cells.  B,  IL-8  KD 
suppresses  cell  invasion  in  both  P  and  TamR  cells.  C,  The  induced  cell 
invasiveness  in  600MPE-P  cells  with  FOXA1-OE  was  reduced  by  simultaneous 
IL-8  KD. 


CHANGES/PROBLEMS 

Nothing  to  Report. 


PRODUCTS 

1.  Publications,  conference  papers,  and  presentations 

Publications 

Fu  X,  Jeselsohn  R,  Pereira  R,  Hollingsworth  EF,  Creighton  CJ,  Li  F,  Shea  M,  Nardone  A,  Angelis  CD, 
Heiser  LM,  Anur  P,  Wang  N,  Grasso  CS,  Spellman  P,  Griffith  OL,  Tsimelzon  A,  Gutierrez  C,  Huang  S, 
Edwards  DP,  Rimawi  MF,  Lopez-Terrada  D,  Hilsenbeck  SG,  Gray  JW,  Brown  M,  Osborne  CK,  Schiff  R. 
FOXA1  overexpression  mediates  endocrine  resistance  by  altering  the  ER  transcriptome  and  IL-8  expression 
in  ER-positive  breast  cancer.  Proc  Natl  Acad  Sci  USA,  Published  online  before  print  October  6,  2016,  doi: 
10.1073/pnas. 16128351 13. 

Conference  presentation 

Xiaoyong  Fu,  Resel  Pereira,  Dongyu  Zhao,  Sung  Yun  Jung,  Rinath  Jeselsohn,  Chad  J.  Creighton,  Martin 
Shea,  Agostina  Nardone,  Carmine  De  Angelis,  Anna  Tsimelzon,  Tao  Wang,  Carolina  Gutierrez,  Shixia 
Huang,  Dean  P.  Edwards,  Mothaffar  F.  Rimawi,  Susan  G.  Hilsenbeck,  Myles  Brown,  Kaifu  Chen,  C.  Kent 
Osborne,  Rachel  Schiff.  FOXA1  induces  a  pro-metastatic  secretome  through  ER-dependent  and  independent 
transcriptional  reprogramming  in  endocrine-resistant  breast  cancer.  Poster  discussion,  San  Antonio  Breast 
Cancer  Symposium,  December  6-10,  2016;  San  Antonio,  TX 

2.  Websites  or  other  Internet  sites 

Nothing  to  Report. 

3.  Technologies  or  techniques 

Nothing  to  Report. 

4.  Inventions,  patent  applications,  and/or  licenses 
Nothing  to  Report. 

5.  Other  products 

a,  The  MCF7L-P  and  TamR  cells  exome  sequencing  data  have  been  deposited  in  the  NCBI  Sequence  Read 
Archive  (SRA)  database  (accession  no.  SRP066629). 

b,  The  related  RNA  sequencing  and  FOXA1  ChIP  sequencing  data  have  been  deposited  in  the  NCBI  Gene 
Expression  Omnibus  (GEO)  database  (SuperSeries  accession  no.  GSE75372). 

c,  Dox-inducible  FOXA1-OE  and  IL-8-KD  breast  cancer  MCF7L  cell  derivatives 

PARTICIPANTS  &  OTHER  COLLABORATIING  ORGANIZATIONS 
1.  Participants 

Name:  Xiaoyong  Fu 

Project  Role:  Principal  Investigator 

Nearest  person  month  worked:  10 

Contribute  to  Project:  Dr.  Fu  has  overseen  this  project,  coordinated  the  efforts  from  collaborators,  worked  on 
experimental  design,  performed  the  experiments  and  data  analysis,  and  written  the  report  and  manuscript. 

Name:  Resel  Pereira 

Project  Role:  Research  Assistant 

Nearest  person  month  worked:  12 

Contribute  to  Project:  Ms.  Pereira  has  assisted  Dr.  Fu  to  perform  the  experiments  and  in  charge  of 
maintaining  our  preclinical  Endo-R  cell  models. 

Name:  Susan  Hilsenbeck 
Project  Role:  Key  personnel 
Nearest  person  month  worked:  1 


8 


Contribute  to  Project:  Dr.  Hilsenbeck  has  participated  in  experimental  design  and  performed  the  biostatistics 
analysis  of  experimental  results. 

Name:  Chad  Creighton 
Project  Role:  Key  personnel 
Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Creighton  has  performed  bioinformatics  analysis  of  next-generation  sequencing 
data. 

Name:  Rinath  Jeselsohn 
Project  Role:  Key  personnel 
Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Jeselsohn  has  collaborated  in  the  RNA-Seq  and  ChIP-Seq  experiments  and  data 
analysis. 

Name:  Rachel  Schiff 
Project  Role:  Key  personnel 
Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Schiff  has  collaborated  in  experimental  design  and  data  interpretation. 

Name:  Tao  Wang 

Project  Role:  Non-Key  personnel 

Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Wang  has  assisted  in  animal  experimental  design  and  data  analysis. 

Name:  Carolina  Gutierrez 
Project  Role:  Non-Key  personnel 
Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Gutierrez  has  assisted  in  development  of  FOXA1/IL-8  IHC  and  data  analysis. 

Name:  Emporia  Hollingsworth 
Project  Role:  Collaborator 
Nearest  person  month  worked:  2 

Contribute  to  Project:  Ms.  Hollingsworth  has  performed  the  FOXA1  FISH  assay  and  data  analysis. 

Name:  Dolores  Lopez-Terrada 
Project  Role:  Collaborator 
Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Lopez-Terrada  has  assisted  in  FISH  data  analysis. 

Name:  Agostina  Nardone 
Project  Role:  Collaborator 
Nearest  person  month  worked:  2 

Contribute  to  Project:  Dr.  Nardone  has  assisted  in  animal  experiment. 

Name:  Martin  Shea 
Project  Role:  Collaborator 
Nearest  person  month  worked:  3 

Contribute  to  Project:  Dr.  Shea  has  assisted  in  animal  experiment. 

Name:  Fugen  Li 

Project  Role:  Collaborator 

Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Li  has  performed  the  bioinformatics  analysis  of  RNA-Seq  and  ChIP-Seq  data. 

Name:  Laura  Heiser 
Project  Role:  Collaborator 
Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Heiser  has  performed  in  the  Exome-Seq  of  Endo-R  cell  models. 
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Name:  Pavana  Anur 
Project  Role:  Collaborator 
Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Anur  has  assisted  in  the  analysis  of  Exome-Seq  data. 

Name:  Nicholas  Wang 
Project  Role:  Collaborator 
Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Wang  has  assisted  in  the  Exome-Seq  data  analysis. 

Name:  Catherine  S.  Grasso 
Project  Role:  Collaborator 
Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Grasso  has  assisted  in  the  next-generation  sequencing  data  analysis. 

Name:  Paul  Spellman 
Project  Role:  Collaborator 
Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Spellman  has  assisted  in  the  Exome-Seq  data  analysis. 

Name:  Joe  Gray 

Project  Role:  Collaborator 

Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Gray  has  assisted  in  the  Exome-Seq  data  analysis. 

Name:  Myles  Brown 
Project  Role:  Collaborator 
Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Brown  has  assisted  in  ChIP-Seq  experimental  design  and  data  analysis. 

Name:  C.  Kent  Osborne 
Project  Role:  Collaborator 
Nearest  person  month  worked:  1 

Contribute  to  Project:  Dr.  Osborne  has  participated  in  experimental  design  and  data  analysis. 

2.  There  are  no  changes  in  the  active  other  support  of  the  PI  or  key  personnel. 

3.  Other  organizations 

Organization  Name:  Texas  Children’s  Hospital 
Location  of  Organization:  Houston,  Texas 

Partner’s  contribution  to  the  project:  Collaboration  in  the  FOXA1  FISH  assay  and  data  analysis 

Organization  Name:  Dana-Farber  Cancer  Institute 
Location  of  Organization:  Boston,  Massachusetts 

Partner’s  contribution  to  the  project:  Collaboration  in  the  FOXA1  ChIP-Seq  and  data  analysis 

SPECIAL  REPORTING  REQUIREMENTS 

None 

APPENDICES 

1 .  Copy  of  the  PNAS  paper 

2.  Abstract  of  San  Antonio  Breast  Cancer  Symposium,  Poster  Discussion,  Dec,  2016 
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Forkhead  box  protein  A1  (FOXA1)  is  a  pioneer  factor  of  estrogen 
receptor  a  (ER)-chromatin  binding  and  function,  yet  its  aberration 
in  endocrine-resistant  (Endo-R)  breast  cancer  is  unknown.  Here, 
we  report  preclinical  evidence  for  a  role  of  FOXA1  in  Endo-R  breast 
cancer  as  well  as  evidence  for  its  clinical  significance.  FOXA1  is 
gene-amplified  and/or  overexpressed  in  Endo-R  derivatives  of  sev¬ 
eral  breast  cancer  cell  line  models.  Induced  FOXA1  triggers  onco¬ 
genic  gene  signatures  and  proteomic  profiles  highly  associated 
with  endocrine  resistance.  Integrated  omics  data  reveal  IL8  as 
one  of  the  most  perturbed  genes  regulated  by  FOXA1  and  ER 
transcriptional  reprogramming  in  Endo-R  cells.  IL-8  knockdown  in¬ 
hibits  tamoxifen-resistant  cell  growth  and  invasion  and  partially 
attenuates  the  effect  of  overexpressed  FOXA1.  Our  study  high¬ 
lights  a  role  of  FOXA1  via  IL-8  signaling  as  a  potential  therapeutic 
target  in  FOXA1 -overexpressing  ER-positive  tumors. 

FOXA1  |  estrogen  receptor  |  breast  cancer  |  transcriptional  reprogramming  | 
endocrine  resistance 

About  75%  of  breast  cancers  express  estrogen  receptor  a 
(ER),  which  is  a  strong  driver  and  therapeutic  target  for 
these  ER-positive  (+)  tumors.  Endocrine  therapy  with  aromatase 
inhibitors  lowers  the  level  of  estrogen;  selective  ER  modulators 
such  as  tamoxifen  (Tam)  bind  to  and  block  ER,  and  down-reg¬ 
ulators  such  as  fulvestrant  (Ful)  bind  to  ER  and  induce  its 
degradation.  Endocrine  therapy  prolongs  disease-free  and 
overall  survival  when  used  in  the  adjuvant  setting  and  can  induce 
long-term  remission  in  some  patients  in  the  metastatic  setting. 
Despite  the  overall  success  of  endocrine  therapy,  tumors  in  more 
than  50%  of  patients  with  metastatic  disease  fail  to  respond,  and 
nearly  all  metastatic  patients  with  initially  responding  tumors 
eventually  experience  tumor  relapse  and  die  from  acquired  re¬ 
sistance  (1,  2).  Although  there  are  many  causes  for  resistance, 
the  most  predominant  mechanisms  include  altered  ER  signaling 
and  interactions  between  ER,  its  coregulators,  and  various 
growth  factor  pathways.  These  alterations  facilitate  adaptation 
from  ligand-dependent  to  ligand-independent  ER  activation, 
which  is  further  triggered  by  cross-talk  with  growth  factor  receptor 
(GFR)  signaling  pathways  (3-6).  However,  the  key  mediators 
of  ER  transcriptional  reprogramming  in  promoting  endocrine- 
resistant  (Endo-R)  breast  cancer  remain  poorly  understood. 

Recently,  a  potential  role  of  the  forkhead  box  protein  A1 
(FOXA1)  has  been  suggested  in  mediating  endocrine  resistance  in 
breast  cancer  (7,  8).  FOXA1  is  termed  a  “pioneer  factor”  because  it 
binds  to  highly  compacted  or  “closed”  chromatin  via  a  domain 
similar  to  that  of  linker  histones  and,  through  its  C-terminal  domain, 


renders  these  genomic  regions  more  accessible  to  other  tran¬ 
scription  factors,  such  as  ER  (9),  progesterone  receptor  (PR) 
(10),  and  androgen  receptor  (AR)  (11).  As  such,  FOXA1  has  a 
key  role  in  demarcating  the  tissue-specific  binding  sites  of  these 
nuclear  receptors  (12).  Together  with  ER,  FOXA1  contributes 
to  the  pattern  of  gene  transcription  that  induces  luminal  cell 
differentiation  (13)  and  represses  the  basal  phenotype  (14).  Like 
ER,  FOXA1  is  associated  with  luminal  subtype  and  good  prog¬ 
nosis  in  breast  cancer  (15,  16).  However,  FOXA1  and  ER  have 
also  been  found  to  be  coexpressed  at  high  levels  in  breast  cancer 
metastases  that  are  resistant  to  endocrine  therapy  (8),  suggesting 
a  continuing  and  potentially  altered  role  of  FOXA1  in  ER+ 
metastatic  and/or  resistant  disease.  A  recent  study  in  endome¬ 
trial  cancer  found  increasing  levels  of  FOXA1  in  metastases, 
even  though  high  levels  of  FOXA1  in  primary  tumors  were 


Significance 

One  of  the  mechanisms  of  endocrine  resistance  in  estrogen  re¬ 
ceptor  a  (ER)-positive  (+)  breast  cancer  is  the  cross-talk  between 
the  ER  and  growth  factor  receptor  pathways  leading  to  altered 
ER  activity  and  a  reprogrammed  ER-dependent  transcriptome. 
However,  key  mediators  of  this  ER-dependent  transcriptional 
reprogramming  remain  elusive.  Here  we  demonstrate  that  fork- 
head  box  protein  A1  (FOXA1)  up-regulation  via  gene  amplifica¬ 
tion  or  overexpression  contributes  to  endocrine  resistance  and 
increased  invasiveness  phenotypes  by  altering  the  ER-dependent 
transcriptome.  We  further  show  that  IL-8,  one  of  the  top  altered 
FOXA1/ER  effectors,  plays  a  key  role  in  mediating  these  pheno¬ 
types  and  is  a  potential  target  to  treat  ER+/FOXA1-high  breast 
cancer.  Our  findings  provoke  a  new  interplay  of  FOXA1  in  the  ER 
transcriptional  program  in  endocrine-resistant  breast  cancer. 
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associated  with  good  outcome  (17).  At  the  molecular  level,  ge¬ 
nome-wide  mapping  of  ds-regulatory  elements  (cistromes)  has 
shown  that  the  FOXAl-binding  motif  is  enriched  in  a  distinct  ER 
cistrome  identified  in  ER+  primary  tumors  from  patients  who  are 
likely  to  relapse,  suggesting  a  functional  link  of  FOXA1  with  ag¬ 
gressive  ER+  disease  (8).  These  contradictory  findings  of  the  sig¬ 
nificance  of  FOXA1  in  early  and  late  tumor  stages  suggest  a 
potentially  dynamic  perturbation  of  FOXA1  in  disease  progression. 
However,  it  remains  unclear  how  FOXA1  is  engaged  in  the  ER 
transcriptional  reprogramming  in  Endo-R  breast  cancer,  and 
whether  there  is  any  aberration  of  FOXA1  that  contributes  to 
this  process. 

The  aim  of  this  study  was  to  evaluate  the  role  of  FOXA1  in 
mediating  endocrine  resistance  in  ER+  breast  cancer  using  a  panel 
of  Endo-R  breast  cancer  cell  line  models,  publicly  deposited 
preclinical  and  clinical  datasets,  and  functional  studies.  Our  hy¬ 
pothesis  was  that  increased  expression  of  FOXA1  in  breast  tumors 
might  contribute  to  endocrine  resistance  and  tumor  progression. 
We  found  that  FOXA1  expression  was  increased  in  several  dif¬ 
ferent  ER+  Endo-R  derivative  cell  lines  compared  with  their  ER+ 
parental  (P)  cells.  Induced  overexpression  of  FOXA1  in  the  P  cells 
elicited  gene  signatures  and  proteomic  profiles  associated  with 
multiple  oncogenic  pathways  as  well  as  endocrine  resistance.  High 
levels  of  FOXA1  mRNA  predicted  poor  outcome  in  patients  with 
ER+  tumors  receiving  Tam.  Integrative  analysis  of  cistromic  and 
RNA  sequencing  (seq)  data  suggested  that  IL-8  serves  as  an  im¬ 
portant  mediator  of  the  FOXA1/ER  transcriptional  reprogram¬ 
ming  to  promote  Endo-R  cell  growth  and  invasion.  We  propose 
that  targeting  IL-8  signaling  is  a  promising  strategy  to  treat  ER+ 
tumors  with  high  levels  of  FOXA1. 

Results 

F0XA1  Gene  Amplification  Is  Associated  with  Tam  Resistance  in  ER+ 
Breast  Cancer  Preclinical  Models.  Five  established  Endo-R  cell 
models  showed  a  stable  phenotype  of  sustained  cell  growth  in  the 
presence  of  estrogen  deprivation  (ED)  or  Tam  (Fig.  SI).  Two 
MCF7  Endo-R  cell  models  were  independently  developed  from 
the  ER+  breast  cancer  MCF7-  L  (18)  and  RN  (19)  lines.  Using 
whole-exome-seq,  we  found  that  the  genomic  region  (14q21.1) 
encompassing  only  the  FOXA1  gene  had  the  highest  focal  am¬ 
plification  ratio  in  Tam-resistant  (TamR)  derivatives  compared 
with  P  cells  in  both  MCF7-L  and  RN  models  [log2  copy  number 
(CN)  ratio  of  3.7  and  3.4  in  Fig.  1 A  and  B  and  Fig.  S2A  and  B , 
respectively].  This  FOXA1  gene  amplification  was  found  only  in 
the  MCF7-L/RN  TamR  but  not  the  ED-resistant  (EDR)  de¬ 
rivative.  Furthermore,  at  a  single  cell  level  there  was  a  highly 
enriched  cell  population  with  FOXA1  amplification  ( FOXA1  vs. 
reference  foci  ratio  >4)  revealed  by  FISH  in  the  MCF7-L/RN 
TamR  compared  with  P  cells  (Fig.  1  C  and  D  and  Fig.  S2  C  and 
D).  Even  in  the  MCF7-L/RN  P  cells,  we  found  a  mixed  cell 
population  with  over  50%  of  cells  showing  a  ratio  >2,  suggesting 
some  level  of  FOXA1  CN  gain  (CNG)  preexisting  in  the  P  cells 
before  developing  endocrine  resistance.  FOXA1  gene  amplifi¬ 
cation  was  also  validated  using  a  genomic  PCR  (gPCR)  assay 
(Fig.  IE).  The  FOXA1-CN  in  MCF7-L/RN  P  cells  was  higher 
than  that  in  the  normal  mammary  epithelial  MCF10A  cells.  In 
fact,  MCF7  cells  had  the  highest  FOXA1-CN  among  a  panel  of 
59  breast  cancer  cell  lines  [data  from  the  Cancer  Cell  Line  En¬ 
cyclopedia  (20)]  (Fig.  S2E).  Two  other  cell  lines,  KPL1  and 
BT474,  also  showed  high  FOXA1-CN.  In  our  gPCR  assay,  we 
also  observed  a  relatively  modest  but  significant  FOXA1-CN 
increase  in  TamR  but  not  EDR  derivatives  of  the  BT474  model. 
FOXA1  amplification  was  not  found  in  two  other  ER+  Endo-R 
models  (ZR75-1  and  600MPE). 

Amplification  of  the  genomic  region  encompassing  the 
FOXA1  gene  has  been  reported  in  primary  and  metastatic  tu¬ 
mors  of  esophagus,  lung,  thyroid,  and  prostate  (21-23).  We  an¬ 
alyzed  the  updated  Cancer  Genome  Atlas  (TCGA)  breast  cancer 
dataset  (n  =  1,105)  for  CN  changes  (24,  25).  Although  FOXA1 
gene  amplification  was  found  only  in  2%  of  all  cases,  20%  of 
tumors  had  FOXA1-CN  alterations  including  both  CNG  and 
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Fig.  1.  FOXA1  gene  amplification  in  preclinical  ER+  Endo-R  cell  models. 
(A)  Overall  CN  across  the  genome  for  MCF7L-TamR  related  to  P  cells.  Log2(CN 
ratio)  is  shown  on  the  vertical  axis.  Each  point  represents  the  log-trans¬ 
formed  CN  ratio  for  each  targeted  exon,  ordered  by  genomic  coordinates 
and  colored  by  chromosome  using  red  and  blue  for  subsequent  chromo¬ 
somes.  Black  lines  show  inferred  segments.  The  arrows  point  to  the  segment 
with  the  highest  focal  amplification  exons.  ( B )  Zoom  into  the  60-Mb  region 
containing  FOXA1  in  A.  CN  ratios  for  the  four  FOXA 1 -targeted  exons  are  shown 
as  purple  squares,  and  the  remaining  targeted  exons  are  red  circles.  The  seg¬ 
ment  containing  FOXA1  contains  just  the  four  targeted  exons  corresponding  to 
FOXA1.  (C)  Representative  images  of  bright-field  and  the  FOXA7-FISH  of 
MCF7L-P  and  TamR  cells  show  the  enrichment  of  gene  amplification  in  MCF7L- 
TamR  vs.  P  cells.  Green  and  red  signals  indicate  the  locations  where  FOXA1  and 
chromosome  14  centromere  reference  (REF)  probes  were  hybridized,  re¬ 
spectively.  [Scale  bar,  100  |im  (bright  field)  and  20  ^im  (FISH).]  (D)  Stacked  bar 
chart  summarizes  the  percentage  of  cells  (n  =  65)  with  FOXAMREF  foci  ratio 
within  indicated  ranges.  ( E )  Normalized  FOXA1-QH  values  from  multiple  Endo-R 
cells  were  calculated  from  the  results  of  a  real-time  gPCR  assay.  The  normal 
diploid  MCFIOAcell  line  was  used  as  the  normalization  control  (CN  =  2,  marked 
by  a  dashed  line).  ( F )  Box-whisker  plots  show  the  FOXA1-CN  across  the  five 
molecular  subtypes  of  breast  cancer  ( n  =  831)  in  the  TCGA  dataset  (24). 


amplification.  The  FOXA1  -CN  was  higher  in  luminal  and  human 
epidermal  growth  factor  receptor  2  (HER2)-enriched  subtypes 
than  in  the  basal  subtype  (Fig.  IF),  which  correlates  with  the 
expression  pattern  of  FOXA1  mRNA  across  the  subtypes  (Fig. 
S3/4).  There  were  more  tumors  in  the  luminal  B  (42%)  with 
FOX47-CNG  and  amplification  than  in  the  luminal  A  (14%) 
subtype,  suggesting  an  association  of  increased  FOXA1-CN  with 
poor  clinical  outcome.  In  partial  support  of  this,  we  found  that 
FOXA1-CN  was  significantly  higher  in  lymph-node  metastases 
compared  with  the  matched  primary  ER+  luminal  tumors  (n  = 
22)  in  a  Gene  Expression  Omnibus  (GEO)  dataset  (accession  no. 
GSE56765)  (26)  (Fig.  S3  B  and  C).  Altogether,  our  preclinical 
data  and  the  clinical  evidence  support  a  hypothesis  that  high 
levels  of  FOXA1-CN  in  aggressive  luminal  tumors  favor  the 
outgrowth  of  Endo-R  tumors  through  a  subclonal  selection  or 
enrichment  in  response  to  endocrine  therapy. 
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F0XA1  Is  Overexpressed  in  Endo-R  Derivatives  and  Is  Essential  for 
Both  P  and  Endo-R  Cell  Growth  in  Multiple  Predinical  Cell  Models. 

Although  FOXA1  amplification/CNG  was  seen  only  in  MCF7-L/RN 
and  BT474  TamR  derivatives,  FOXA1  mRNA  levels  were 
higher  in  the  TamR  derivatives  than  in  the  P  cells  of  all  five 
models  (MCF7-L/RN,  BT474,  ZR75-1,  and  600MPE)  measured 
by  quantitative  reverse-transcription  (qRT)-PCR  (Fig.  24). 
Similarly,  increased  FOXA1  mRNA  was  also  observed  in  the 
EDR  derivatives  of  ZR75-1,  600MPE,  and  BT474  models.  In¬ 
creased  FOXA1  protein  levels  measured  by  Western  blot  were 
observed  in  the  Endo-R  derivatives  compared  with  their  P  cells 
in  all  five  models  (Fig.  2 B).  ER  protein  was  retained  in  all  but 
one  of  the  Endo-R  cell  lines  compared  with  P  cells;  the  ZR75-1 
Endo-R  model  had  no  detectable  ER.  Protein  levels  of  classical 
ER-regulated  genes  such  as  PGR  and  BCL2,  as  well  as  GATA3 
(encoding  GATA-binding  protein  3),  which  also  regulates  ER 
expression  (27),  were  down-regulated  in  most  of  these  Endo-R 
derivatives  compared  with  P  cells  (Fig.  2 B  and  Fig.  S4  ^4-C), 
suggesting  a  continuous  blockade  of  the  classical  ER  transcrip¬ 
tional  program  that  is  also  seen  in  our  previously  reported  Endo-R 
xenograft  mouse  model  (5).  Importantly,  high  FOXA1  protein 
levels  were  also  observed  by  immunohistochemistry  (IHC)  in 
acquired  Endo-R  MCF7L  xenograft  tumors  in  vivo  compared 
with  estrogen-treated  controls  (Fig.  2C). 

To  determine  the  role  of  ER  and  FOXA1  in  endocrine  re¬ 
sistance,  we  evaluated  cell  growth  of  various  P  and  Endo-R  de¬ 
rivatives  in  response  to  two  validated  siRNAs  targeting  ER  and 
FOXA1  (Fig.  2D).  Knocking  down  ER  in  the  MCF7-L/RN 
models  significantly  inhibited  both  P  and  Endo-R  cell  growth 
(Fig.  2  E  and  F).  Both  ZR75-1-P  and  600MPE-P  cells  were  also 
sensitive  to  ER  knockdown;  however,  cell  growth  was  affected  to 
a  lesser  extent  by  ER  knockdown  in  their  Endo-R  derivatives 


(Fig.  2  G  and  H).  However,  FOXA1  knockdown  substantially 
inhibited  the  growth  of  P  and  Endo-R  derivatives  of  all  pre- 
clinical  models,  suggesting  an  important  role  of  FOXA1  on 
breast  cancer  cell  growth  even  in  the  setting  of  endocrine  re¬ 
sistance  and  even  in  resistant  cells  that  are  not  affected  by 
ER  knockdown. 

F0XA1  Overexpression  Elicits  an  Endo-R  Gene  Signature  and  Predicts 
Poor  Outcome  in  Patients  with  ER+  Tumors.  To  better  understand  the 
role  of  increased  FOXA1  in  Endo-R  cells,  we  established  a  stable 
MCF7L/FOXA1  cell  model  with  doxycycline  (Dox)-inducible 
FOXA1  overexpression.  The  extent  of  FOXA1  overexpression  in 
the  MCF7L/FOXA1  cells  after  Dox  induction  vs.  without  Dox  was 
comparable  to  that  observed  in  the  MCF7L-TamR  vs.  P  cells  (Fig. 

3/4).  RNA-seq  analysis  revealed  a  total  of  440  genes  up-regulated 
and  217  genes  down-regulated  [|Gfold|  (28)  >1.5,  false  discovery 
rate  (FDR)  <0.05]  in  -l-Dox  vs.  -Dox  cells  (Fig.  3 B).  Functional 
annotation  of  these  up-regulated  genes  in  the  Database  for  An¬ 
notation,  Visualization  and  Integrated  Discovery  (DAVID)  (29) 
showed  a  robust  enrichment  of  Gene  Ontology  (GO)  terms  that 
included  “cell  motion  and  migration,”  “response  to  hypoxia,”  and 
“blood  vessel  development”  (P  <  0.001).  Interestingly,  within  the 
down-regulated  genes,  the  most  enriched  GO  term  was  “response 
to  estrogen”  ( P  =  0.0015),  suggesting  a  reduction  of  ligand- 
dependent  classic  ER  transcriptional  activity  in  this  model,  which 
could  be  partially  due  to  the  modest  decrease  of  ER  expression  © 
itself  (Fig.  3 B,  Lower).  We  further  used  Gene  Set  Enrichment  § 

Analysis  (GSEA)  (30)  to  interrogate  the  oncogenic  gene  signa-  - 

tures  from  MSigDB  (31).  The  MCF7L/FOXA1  gene  expression  " 
profile  was  highly  correlated  to  the  gene  sets  enriched  in  MCF7 
cells  overexpressing  ligand-activated  epidermal  GFR  or  constitu- 
tively  active  MEK1,  or  in  epithelial  cell  lines  over  expressing  an 
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Fig.  2.  Increased  FOXA1  expression  in  multiple 
Endo-R  cell  models.  (A)  FOXA1  mRNA  levels  de¬ 
termined  by  qRT-PCR  across  various  Endo-R  cell 
models.  The  amount  of  FOXA1  mRNA  from  MCF7L-P 
cells  was  set  as  normalization  control  (=  1). 
(B)  FOXA1,  ER,  and  PR  protein  expression  using  se¬ 
lected  antibodies  in  Western  blot  analysis  across  vari¬ 
ous  Endo-R  cell  models.  GAPDH  was  used  as  a  loading 
control.  (C)  Scatter  dot  plots  of  FOXA1  Allred  score  in 
MCF7L  Endo-R  xenograft  tumors  measured  by  IHC. 
Xenograft  tumors  in  ovariectomized  nude  mice  with 
estrogen  pellets  (E2),  or  without  E2  but  treated  with 
Tam  or  without  Tam  (ED),  were  harvested  when  the 
tumor  volume  reached  1,000  mm3.  Data  are  repre¬ 
sented  as  mean  ±  SEM  (n  >  5  in  each  group).  ***P  < 
0.001,  Bonferroni  post  hoc  test  (multiple  testing- 
corrected).  (D)  Western  blots  showing  the  protein 
levels  of  FOXA1  and  ER  in  MCF7L-TamR  cells  with 
gene  knockdown.  ( E-H )  Cell  growth  within  a 
6-d  period  after  siRNA  knockdown  of  nonspecific 
(N.S.)  target,  ER,  or  FOXA1  in  MCF7L,  MCF7RN, 
ZR75-1 ,  and  600MPE  Endo-R  cell  models.  Cell 
growth  in  N.S.  knockdown  was  used  as  normaliza¬ 
tion  control  (100%).  Data  represent  means  ±  SEM, 
*P  <  0.05,  **P  <  0.01,  ***P  <  0.001,  two-sided  t  test 
for  all  comparisons  between  N.S.  and  ER/FOXA1 
knockdown. 
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oncogenic  KRAS,  suggesting  that  FOXA1  overexpression  en¬ 
hances  GFR  downstream  signaling.  In  addition,  this  FOXA1- 
induced  transcriptomic  profile  was  significantly  enriched  for  the 
gene  set  that  was  up-regulated  in  the  MCF7  xenograft  tumors  that 
acquired  resistance  to  multiple  endocrine  therapies  from  our 
previously  published  study  (3)  (Fig.  3C).  These  data  suggest  that 
increased  FOXA1  potentially  drives  a  transcriptional  program 
associated  with  high  GFR  signaling  that  contributes  to  tumor 
aggressiveness  and  endocrine  resistance. 

Because  differentially  expressed  genes  from  our  FOXAl-over- 
expressing  MCF7L/FOXA1  preclinical  cell  model  were  enriched 
for  genes  in  our  previously  described  signature  from  Endo-R  xe¬ 
nograft  models  (3),  we  asked  whether  FOXA1  levels  were  corre¬ 
lated  with  the  endocrine  resistance  signature  score  in  clinical 
samples.  Indeed,  high  FOXA1  mRNA  levels  in  752  ER+  tumors 
(32)  were  positively  correlated  with  the  Endo-R  gene  signature 
(Fig.  3D,  Spearman  correlation,  r  =  0.083,  P  =  0.011).  Next,  we 
tested  the  endocrine  response  in  our  Dox-induced  FOXA1 -over- 
expressing  MCF7L  and  ZR75-1  cell  models.  The  highest  levels  of 
FOXA1  induced  endocrine  resistance  in  both  cell  models  (Fig.  3  E 
and  F).  Specifically,  increased  FOXA1  expression  was  significantly 
associated  with  decreased  endocrine  sensitivity  to  Tam  in  MCF7L/ 
FOXA1  cells  and  to  ED  in  ZR75-1/FOXA1  cells,  in  a  FOXA1 
level-dependent  manner.  The  role  of  FOXA1  expression  levels  in 
treatment  response  was  also  reflected  in  clinical  samples.  In  a 
metaanalysis  of  published  datasets  (kmplot.com)  (33),  we  found 
that  the  top  quartile  of  FOXA1  mRNA  levels  was  associated  with 
poor  relapse-free  survival  (RFS)  in  patients  with  ER+  tumors  re¬ 
ceiving  Tam  ( n  =  615,  P  =  0.029),  but  not  in  patients  without  en¬ 
docrine  therapy  (n  =  500,  P  =  0.81)  (Fig.  3  G  and  H).  Collectively, 
these  data  suggest  that  high  FOXA1  expression  is  functionally, 
biologically,  and  clinically  associated  with  endocrine  resistance. 

Proteomic  Profiles  Perturbed  by  F0XA1  Overexpression  Are  Associated 
with  Multiple  Oncogenic  Pathways.  Because  of  the  clinical  evidence 
for  the  potential  role  of  FOXA1  in  mediating  endocrine  resistance, 
we  wanted  to  further  dissect  its  downstream  signaling  pathways. 
For  this,  we  applied  reverse-phase  protein  arrays  (RPPA)  to  de¬ 
termine  the  proteomic  changes  in  our  FOXAl-overexpressing  ER+ 
cell  models,  using  a  total  of  204  validated  antibodies.  Proteins 
differentially  expressed  between  -i-Dox  (at  day  2  or  5)  and  -Dox 
samples  were  identified  (Dataset  SI,  one-way  ANOVA,  P  <  0.05) 
and  visualized  in  heat  maps  following  hierarchical  clustering  (Fig.  4 
A-C).  Consistent  with  the  RNA-seq  data,  the  protein  levels  of  ER 
and  the  products  of  its  classically  regulated  genes  (e.g.,  PGR, 
BCL2,  and  MYC)  were  decreased  in  the  MCF7L/FOXA1  +D0X 
cells  (Fig.  S4D).  Assigning  the  total  proteins  assessed  by  RPPA  into 
Kyoto  Encyclopedia  of  Genes  and  Genomes  (KEGG)  (34)  -de¬ 
fined  cancer  pathways,  we  tracked  the  pathway  activation  status  by 
comparing  the  averaged  signals  within  each  pathway  between 
-/+  Dox  samples.  We  found  that  the  GFR  pathways  of  focal  adhe¬ 
sion,  ERBB2,  and  insulin  were  overactivated  in  both  the  MCF7L/ 
FOXA1  +Dox  and  ZR75-1/FOXA1  +Dox  cells  (Fig.  4  D  and  E, 
P  <  0.001).  The  NOTCH  pathway,  which  previously  has  been 
shown  to  be  overactivated  in  Endo-R  breast  cancer  cells  (35),  did 
not  seem  significantly  perturbed  by  FOXA1  overexpression  in  our 
cell  models,  possibly  due  to  the  low  number  of  representative 
pathway  proteins  in  this  RPPA  assay.  The  decreased  ER  and  in¬ 
creased  GFR  downstream  signaling  in  the  MCF7L/FOXA1  -l-Dox 
cells  was  further  confirmed  by  Western  blot  showing  a  FOXA1- 
dependent  effect  (Fig.  S5A).  The  600MPE/FOXA1  +D0X  cells 
showed  less  enhanced  GFR  signaling,  possibly  due  to  an  endoge¬ 
nously  hyperactivated  mitogen-activated  protein  kinase  (MAPK) 
pathway  caused  by  a  KRAS  mutation  in  this  line  (36)  (Fig.  4 F). 
Overall,  there  were  23  commonly  up-regulated  and  1  down-regulated 
(GATA3)  proteins  across  all  three  cell  models  upon  FOXA1  over¬ 
expression  (Fig.  4  G  and  H).  The  significantly  enriched  insulin  and 
mechanistic  target  of  rapamycin  (mTOR)  pathways  represented  by 
the  23  commonly  up-regulated  proteins  (Fig.  41),  together  with  the 
commonly  decreased  luminal  lineage  marker  and  reciprocal 
ER  regulator  GAT  A3  (27),  further  support  the  role  of  increased 
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Fig.  3.  FOXA1  overexpression  is  associated  with  endocrine  resistance  and 
poor  clinical  outcome.  (A)  FOXA1  IHC  of  MCF7L-P  and  TamR,  and  MCF7L7 
FOXA1  ±  Dox  cells.  (Scale  bar,  100  |am.)  ( B )  Heat  map  of  differentially  expressed 
genes  (|Gfold|  >1.5)  after  FOXA1  overexpression  in  MCF7L/FOXA1  cells.  The 
levels  of  ER  protein  in  ±  Dox  cells  are  shown  at  the  bottom  panel.  DAVID 
functional  annotation  shows  the  enriched  GO  terms  for  the  genes  up-  or 
down-regulated  upon  FOXA1  overexpression.  (C)  GSEA  shows  the  correlation 
of  the  gene  expression  profile  of  MCF7L/FOXA1  cells  with  a  gene  set  derived 
from  the  Endo-R  xenograft  tumors  (3).  (D)  Spearman  correlation  of  FOXA1 
mRNA  levels  and  the  signature  scores  of  the  Endo-R  gene  set  in  ER+  breast 
tumors  (n  =  752).  (£  and  F)  Cell  growth  was  measured  in  MCF7L/FOXA1  and 
ZR75-1/FOXA1  cells  with  0,  0.25,  or  1  |ig/mL  Dox  in  presence  of  E2  or  various 
endocrine  therapies.  E2-treated  cells  were  used  as  normalization  controls  for 
anti-estrogen  groups  (ED,  Tam,  and  Ful).  The  induced  FOXA1  proteins  are 
shown  by  Western  blots.  Data  represent  means  ±  SEM,  *P  <  0.05,  **P  <  0.01, 
two-sided  t  test  for  all  indicated  comparisons.  (G)  Kaplan-Meier  plots  show  RFS 
in  ER+  patients  receiving  Tam  but  without  chemotherapy  (n  =  615),  who  were 
stratified  by  FOXA1  mRNA  levels  at  the  top  quartile  (25%)  vs.  the  rest  (75%). 
(H)  Same  analysis  in  ER+  patients  without  endocrine  treatment  (n  =  500).  P 
value  was  calculated  by  using  the  log-rank  test.  Analysis  was  performed  using 
an  online  tool  and  resource  at  kmplot.com/analysis  (33). 


FOXA1  in  augmenting  GFRs  and  suppressing  the  classical  ER 
signaling  in  ER+  breast  cancer. 

We  also  performed  RPPA  analysis  in  the  MCF7L-TamR  cells 
with  FOXA1  knockdown.  Interestingly,  the  level  of  proteins  re¬ 
lated  to  the  classical  ER  pathway  such  as  PR  and  GAT  A3,  which 
was  decreased  in  TamR  vs.  P  cells,  was  restored  by  FOXA1 
knockdown  (Fig.  S 5B).  Furthermore,  FOXA1  knockdown  in 
MCF7L-TamR  cells  suppressed  the  oncogenic  pathways  (e.g., 
ERBB2  and  insulin  receptor)  that  otherwise  were  enhanced  in 
FOXAl-overexpressing  P  cells  (Fig.  S5C).  The  overall  proteomic 
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changes  in  the  P  cells  with  FOXA1  overexpression  were  inversely 
correlated  to  the  changes  in  the  TamR  cells  with  FOXA1 
knockdown  (Fig.  S 5D,  Pearson  correlation,  r  =  -0.645,  P  = 
0.017).  Together  with  previous  transcriptomic  data,  these  find¬ 
ings  point  to  a  dominant  role  of  increased  FOXA1  in  aug¬ 
menting  oncogenic  signaling  pathways  in  endocrine  resistance, 
resulting  in  an  inhibitory  effect  on  ER  expression  and  classic  ER 
transcriptional  activity. 

An  Integrative  Approach  Identifies  IL8  As  One  of  the  Most  Perturbed 
Genes  Regulated  by  FOXA1  in  Endo-R  Cells.  To  further  investigate 
the  direct  impact  of  FOXA1  on  gene  expression,  we  performed 
FOXA1  genome-wide  chromatin  immunoprecipitation  followed 
by  high-throughput  seq  (ChIP-seq)  in  MCF7L-P  and  TamR  cells. 
A  total  of  37,227  and  53,215  FOXA1  binding  events  were  found 
in  MCF7L-P  and  TamR  cells,  respectively  (Fig.  S6).  Among 
these  binding  events,  there  were  21,449  shared  FOXA1  binding 
events,  which  accounted  for  58%  and  40%  of  total  binding 
events  in  P  and  TamR  cells,  respectively.  Within  the  distinct 
binding  events  in  P  and  TamR  cells,  the  highest  enrichment  was 
the  FOXA1  motif,  followed  by  the  GATA  motif  in  P  cells,  and 
the  BCL11A  and  JUN/FOS  motifs  in  TamR  cells,  suggesting 


significant  FOXAl-chromatin  binding  in  both  P  and  TamR  cells, 
albeit  on  different  sites.  In  parallel  with  the  cistromic  profiling, 
we  also  obtained  the  transcriptomic  profiles  of  both  MCF7L-P 
and  TamR  cells  using  RNA-seq.  In  an  effort  to  identify  the 
downstream  signaling  associated  with  FOXA1  in  endocrine  re¬ 
sistance,  we  integrated  the  RNA-seq  transcriptomic  data  with 
the  FOXA1  ChIP-seq  data  described  above.  The  genes  preferen¬ 
tially  expressed  in  either  TamR  or  P  cells  tended  to  have  more 
FOXA1  binding  events  (tags)  represented  by  reads  per  million  per 
nucleotide  (RPM),  supporting  the  notion  that  FOXA1  is  indeed 
important  for  defining  the  distinct  gene  patterns  in  both  TamR 
and  P  cells  (Fig.  5 A). 

Next,  we  focused  on  the  top  genes  that  are  highly  expressed  in 
MCF7L-TamR  compared  with  P  cells  and  that  also  carry  the 
most  abundant  FOXA1  tags  (RPM)  in  TamR  cells  around  their 
gene  regions  (log2  ratio  >0.5)  (Fig.  5B).  The  enriched  GO  terms 
within  these  top  genes  include  “blood  vessel  development”  ( IL8 , 
CTGF ,  FOX ,  ROBOl,  HEY1 ,  and  GBX2 )  and  “cell  migration” 
(IL8,  CTGF ,  ROBOl ,  GBX2 ,  and  NR2F1 ),  reminiscent  of  the 
GO  terms  enriched  in  the  FOXAl-overexpressing  MCF7L-P 
cells.  Indeed,  50%  of  these  genes  (e.g.,  IF8,  CTGF ,  and  FOX) 
were  highly  up-regulated  in  MCF7L/FOXA1  -l-Dox  vs.  -Dox 


ZR75-1/FOXA1 


A  A 
vs?  <? 


Fig.  4.  FOXA1  overexpression  in  ER+  breast  cancer  cell  lines  induces  proteomic  perturbations  in  multiple  oncogenic  signaling  pathways.  Heat  maps  of  RPPA  data 
representing  differentially  expressed  proteins  (one-way  ANOVA,  P  <  0.05)  in  MCF7L/FOXA1  (A),  ZR75-1/FOXA1  ( B ),  and  600MPE/FOXA1  (C)  cells  upon  Dox 
addition  for  2  or  5  d.  ( D-F )  Signaling  perturbations  by  FOXA1  overexpression  in  14  cancer-related  KEGG  pathways  were  evaluated  in  these  cell  models  of  A,  B,  and 
C,  by  averaging  the  expression  levels  of  proteins  available  from  RPPA  (numbers  in  parentheses)  within  the  same  pathway,  followed  by  subtraction  of  basal  levels 
in  -Dox  cells.  A  paired  one-sided  t  test  was  applied  and  the  P  value  was  plotted  as  minus  logi0-transformed.  The  perturbations  in  pathways  with  P<  0.05  (P  =  0.05 
is  marked  by  a  dashed  line)  were  statistically  significant.  Venn  diagrams  showing  the  overlapping  proteins  up-  (G)  or  down-regulated  (H)  across  all  three  cell 
models  upon  FOXA1  overexpression  (day  5  of  +DOX  vs.  -Dox,  P  <  0.05).  (/)  Enriched  KEGG  pathways  represented  by  the  commonly  up-regulated  proteins  in  all 
three  cell  models.  Significance  in  enrichment  was  calculated  by  Fisher's  exact  test  (P  =  0.05  is  marked  by  a  dashed  line). 
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cells  (log2  ratio  >0.5),  suggesting  FOXAl-dependent  regulation. 
Furthermore,  we  found  that  there  was  a  significant  overlap  be¬ 
tween  the  genes  highly  represented  in  MCF7L/FOXA1  -l-Dox 
cells  (n  =  440,  Gfold  >1.5,  FDR  <0.05)  and  the  TamR  signature 
genes  (n  =  428,  |log2-ratio|  >1.5,  FDR  <0.05)  (Fig.  5C,  Fisher’s 
exact  test,  P  <  0.0001).  These  genes  included  IL8,  CTGF ,  and 
LOX ,  further  suggesting  that  they  may  play  a  role  in  the  FOXAl- 
dependent  mechanism  of  Tam  resistance.  Finally,  we  found  that 
about  50%  of  genes  (including  IL8 ,  CTGF ,  and  GBX2)  highly 
expressed  in  TamR  cells  with  enhanced  FOXA1  binding  sites 
were  repressed  by  ER  knockdown  in  TamR  cells,  suggesting  that 
at  least  some  of  the  genes  regulated  by  increased  FOXA1  are 
also  dependent  on  ER. 

We  verified  the  robust  increase  in  mRNA  levels  of  IL8,  the 
gene  at  the  top  of  the  list,  in  our  two  independent  TamR  cell 
models  from  the  MCF7  line  (L  and  RN)  by  qRT-PCR  (Fig.  5D). 
In  addition,  significantly  increased  IL-8  expression  was  also  found 
in  both  600MPE  and  ZR75-1  Endo-R  cell  derivatives  compared 
with  their  P  cells,  although  the  magnitude  was  much  smaller  in 
ZR75-1  Endo-R  cells.  Because  ER  expression  is  maintained  in 
these  Endo-R  cells,  except  the  ZR75-1  model,  we  postulated  that 
the  robust  up-regulation  of  IL-8  might  need  both  ER  and  FOXA1. 
It  has  been  reported  that  the  FOXAl-mediated  reprogramming 
of  ER  binding  is  associated  with  the  differential  ER-binding 
program  in  ER+  tumors  from  patients  with  poor  outcome  (8).  We 
hypothesized  that  increased  FOXA1  may  contribute  to  ER 
transcriptional  reprogramming  in  our  Endo-R  cells.  To  better 
appreciate  the  impact  of  increased  FOXA1  on  transcriptional 
switching  of  ER  from  a  ligand-dependent  to  a  growth  factor- 
induced  and  ligand-independent  program,  we  further  integrated  our 
RNA-seq  data  of  FOXAl-overexpressing  MCF7L/FOXA1  cells 
with  the  existing  FOXA1  cistrome  (7)  as  well  as  the  ER  cistrome 
induced  by  estrogen  (E2)  or  epidermal  growth  factor  (EGF)  (37) 
in  MCF7  cells.  As  shown  in  Fig.  5E,  there  were  overlapping  as  well 
as  distinct  subsets  among  the  genes  predicted  from  the  FOXA1 
and  ER  cistromes  [genes  with  FOXA1/ER  binding  sites  ±20  kb  of 
their  transcription  start  sites  (TSS)].  We  intersected  the  up- 
regulated  (UP)  or  not-altered  (NA)  genes  induced  by  FOXA1  in 
our  MCF7L/FOXA1  cells  and  the  list  of  genes  putatively  associ¬ 
ated  with  the  FOXA1  cistrome  and  the  ER  cistrome  induced  by 
EGF,  E2,  or  both.  We  found  that  the  FOXA1-UP  genes  were 
highly  enriched  for  the  genes  associated  with  FOXA1  binding 
and  ER  binding  induced  by  EGF  but  not  by  E2  (Fig.  5F,  Fisher’s 
exact  test,  P  <  0.001).  Notably,  IL8  and  CTGF  were  found  again 
among  the  genes  with  both  FOXA1  binding  and  EGF-induced  ER 
binding,  suggesting  that  the  gene  regulation  by  increased  FOXA1 
involves  a  growth  factor-stimulated  ER-dependent  process.  These 
data  suggest  that  high  levels  of  FOXA1  may  coordinate  in  the  ER 
transcriptional  reprogramming  toward  a  more  growth  factor-induced 
cistromic  profile,  leading  to  endocrine  resistance  by  a  mechanism 
similar  to  that  we  had  previously  reported  in  an  ER/HER2-positive 
MCF7  cell  model  (38,  39). 

Increased  FOXA1,  Together  with  ERf  Coregulates  IL-8  Expression.  Next, 
we  investigated  the  regulation  of  IL-8.  Previous  cistromic  data  in 
MCF7  cells  (7,  37)  revealed  that  there  were  two  FOXA1  binding 
sites  at  the  distal  (dis. )  and  proximal  (pro.)  regions  upstream  of  the 
IL8  TSS,  and  one  EGF-stimulated  ER  binding  site  at  the  dis.  re¬ 
gion  (Fig.  6A).  Our  FOXA1  ChIP-seq  data  showed  that  the 
FOXA1  binding  at  the  dis.  region  of  IL8  in  MCF7L-TamR  cells 
was  enhanced  compared  with  that  in  MCF7L-P  cells  upon  either 
Tam  or  E2  treatment  (Fig.  6 B).  Using  ChIP  followed  by  qPCR,  we 
verified  the  enhancement  of  FOXA1  binding  at  the  dis.  region  in 
MCF7L-TamR  cells  (Fig.  6C).  Furthermore,  there  was  an  en¬ 
hanced  recruitment  of  ER  at  the  dis.  region  in  MCF7L-TamR 
cells,  the  same  region  where  ER  binding  was  previously  shown  in 
MCF7  cells  treated  by  EGF  (37)  or  a  mitogenic  mixture  (8)  in  the 
absence  of  the  E2  ligand.  The  ER  binding  at  the  pro.  region  was 
also  enhanced  in  TamR  vs.  P  cells  in  this  ChIP-qPCR  assay.  These 
data  suggest  that  ER  can  regulate  IL-8  in  a  ligand-independent 
manner  in  the  context  of  high  GFR  signaling  associated  with  en- 
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Fig.  5.  Integrative  analysis  revealed  IL8  as  a  target  of  increased  FOXA1  in  ER 
transcriptional  reprogramming.  (/A)  Integrated  RNA-seq  and  FOXA1  ChIP-seq 
data  in  MCF7L-P  and  TamR  cells.  Genes  aligned  in  RNA-seq  were  calculated  for 
their  expression  log2  ratio  of  fragments  per  kilobase  of  transcript  per  million 
mapped  reads  (FPKM)  in  TamR  vs.  P  cells,  by  which  the  genes  were  sorted  in  a 
descending  order.  FOXA1  binding  events  (tags)  within  ±  20  kb  of  each  gene's  TSS 
were  counted  and  represented  by  average  normalized  RPM  for  every  300  con¬ 
secutive  genes  along  the  order  of  sorted  genes  from  RNA-seq.  These  FOXA1  tags 
were  plotted  separately  for  P  (in  blue)  and  TamR  (in  red)  cells.  ( B )  Heat  maps  of 
genes  with  high  expression  [log2(FPKM)]  and  with  enriched  FOXA1  binding  (log2 
ratio  of  RPM)  in  TamR  vs.  P  cells.  Heat  maps  of  the  expression  in  these  genes  (log2 
ratio  of  FPKM)  in  MCF7L/FOXA1  +Dox  vs.  -Dox  cells  and  in  TamR  cells  with  si-ER 
vs.  si-N.S.  knockdown  are  also  shown.  (C)  Venn  diagram  showing  the  overlap 
genes,  including  IL8,  CTGF,  and  LOX,  between  the  FOXA1 -overexpression  (O.E.) 
up-regulated  genes  and  the  MCF7L-TamR  signature  genes.  P  value  was  calcu¬ 
lated  by  Fisher's  exact  test.  (D)  IL-8  gene  expression  measured  by  qRT-PCR  in  four 
Endo-R  cell  models.  Data  represent  means  ±  SEM,  *P  <  0.05,  **P  <  0.01,  ***P  < 
0.001,  two-sided  t  test  for  all  comparisons  between  Endo-R  and  P  cells.  (£)  Venn 
diagram  showing  the  overlap  among  the  predicted  genes  with  the  binding  of 
FOXA1  and  EGF/E2-stimulated  ER  within  ±  20  kb  of  TSS  in  MCF7  cells.  The  gene 
sets  with  highlighted  numbers  (in  red)  were  used  for  the  following  analysis. 
(E)  The  genes  induced  (UP)  or  not  altered  (NA)  by  FOXA1  overexpression  were 
intersected  with  the  FOXA1  and  ER  cistromes.  Gene  enrichment  within  the 
FOXA1  -UP  gene  set  for  the  genes  associated  with  FOXA1  binding  and  ER  binding 
induced  only  by  EGF,  E2,  or  both  was  compared  with  the  enrichment  within  the 
FOXA1-NA  gene  set.  ***P  <  0.001,  Fisher's  exact  test. 


docrine  resistance.  IL-8  mRNA  levels  in  the  two  TamR  cell  models 
(MCF7L  and  600MPE)  were  reduced  by  either  ER  or  FOXA1 
knockdown  (Fig.  6  D  and  E),  with  the  strongest  reduction  in  the 
MCF7L-TamR  cells  from  knockdown  of  FOXA1,  suggesting  that 
these  binding  events  are  also  biologically  relevant.  In  parallel,  se¬ 
cretory  IL-8  protein  was  dramatically  induced  by  FOXA1  over¬ 
expression  in  MCF7L-P  cells;  the  increased  IL-8  by  FOXA1 
was  substantially  reduced  by  simultaneous  ER  knockdown 
(Fig.  6F).  This  phenomenon  could  be  recapitulated  in  a  second 
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600MPE/FOXA1  cell  model  (Fig.  6 G),  supporting  the  notion 
that  high  levels  of  FOXA1  and  ER  might  coregulate  IF-8  ex¬ 
pression.  In  line  with  the  increased  FOXA1  and  ER  protein 
levels  in  our  MCF7L  Endo-R  xenograft  tumors  (Fig.  2 C  and  Fig. 
S7  A  and  B),  IF-8  expression  was  also  up-regulated  in  both  TamR 
and  EDR  tumors  in  this  xenograft  model  (Fig.  6  H  and  I). 
Moreover,  we  measured  FOXA1  and  IF-8  protein  levels  in  a 
tissue  microarray  composed  of  primary  ER+  tumor  specimens 
(n  =  85)  archived  in  our  tumor  bank.  FOXA1  staining  was  localized 
mainly  in  nuclei,  whereas  IF-8  staining  was  mainly  in  paranuclear 
regions  of  cancer  cells  (Fig.  6J).  Paranuclear  staining  of  IF-8  was 
also  seen  both  in  MCF7F-TamR  xenograft  tumors  and  in  TamR 
cells  in  vitro  (Fig.  S7C).  The  proportion  of  IF-8  positive  tumors 
gradually  increased  with  increasing  FOXA1  in  these  ER+  human 
primary  breast  tumors  (Fig.  6 K,  Fisher’s  exact  test,  P  =  0.006), 
consistent  with  the  notion  that  FOXA1  regulates  IF-8. 

IL-8  Mediates  the  Effect  of  Augmented  FOXA1  on  Cell  Growth, 
Invasion,  and  Endocrine  Resistance.  To  evaluate  the  role  of  IF-8 
as  a  downstream  effector  of  increased  FOXA1  in  endocrine  re¬ 


sistance,  we  obtained  the  transcriptomic  profiles  of  MCF7F- 
TamR  cells  with  FOXA1  or  IF-8  knockdown  by  using  RNA-seq. 
Comparing  the  genes  differentially  expressed  in  TamR  cells 
upon  FOXA1  knockdown  ( |  Gfold  |  >0.5),  we  found  a  striking 
similarity  in  the  expression  patterns  of  these  same  genes  in 
TamR  cells  upon  IF-8  knockdown  (Fig.  S&4),  suggesting  a  cru¬ 
cial  role  of  IF-8  in  the  gene  expression  perturbations  induced 
by  FOXA1  in  Endo-R  cells.  Fike  FOXA1  knockdown,  IF-8 
knockdown  potently  inhibited  cell  growth  in  both  MCF7F-P  and 
TamR  cells,  with  a  greater  growth  inhibitory  effect  in  TamR  than 
in  P  cells  (Fig.  1A,  two-way  ANOVA  interaction  test,  P  <  0.001). 
This  knockdown  effect  was  rescued  by  coexpression  of  an  IF-8 
cDNA  with  the  IF-8  siRNA  sequences  (#2)  targeting  the 
3-UTR  region  of  the  IL8  gene,  but  not  with  the  other  siRNA 
(#1)  targeting  the  coding  region  of  IL-8  (Fig.  IB).  Because 
FOXA1  knockdown  in  MCF7F-TamR  cells  suppressed  multi¬ 
ple  oncogenic  pathways  that  otherwise  were  up-regulated  in 
FOXA1 -over expressing  P  cells  (Fig.  S5C),  we  asked  whether  IF-8 
knockdown  in  TamR  cells  leads  to  a  similar  change  in  signaling. 
Indeed,  the  activated  signaling  of  multiple  GFR  downstream 
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Fig.  6.  Increased  FOXA1  and  ER  regulate  IL-8  expression  in  ER+  breast  cancer.  (A)  Schematic  diagram  of  ER  and  FOXA1  binding  within  the  IL8  gene  locus  as  defined  by 
EGF-stimulated  ER  ChIP-on-chip  (37)  and  FOXA1  ChIP-sEq.  (7)  in  MCF7  cells.  ( B )  Snapshot  of  FOXA1  continuous  peaks  from  ChIP-seq  data  showing  the  binding  pattern 
upstream  of  the  IL8  gene  TSS  in  MCF7L-P  and  MCF7L-TamR  cells  treated  with  Tam  or  E2.  (C)  FOXAI-ChIP  (Left)  and  ER-ChIP  (Right)  followed  by  qPCR  of  binding  regions  in 
MCF7L-P  and  TamR  cells.  Quantification  of  amplified  binding  regions  was  calculated  as  fold  enrichment  by  normalizing  to  an  intergenic  sequence  as  a  negative  control. 
(D  and  E)  Measurement  of  IL-8  mRNA  by  qRT-PCR  in  MCF7L-TamR  and  600MPE-TamR  cells  with  either  ER  or  FOXA1  knockdown.  N.S.,  nonspecific;  #1  and  #2,  two  different 
siRNA  sequences.  (F  and  G)  ELISA  of  IL-8  protein  in  culture  media  of  MCF7L/FOXA1  and  600MPE/FOXA1  cells  -/+Dox  in  the  absence/presence  of  ER  knockdown. 
(H)  Representative  H&E  staining  and  IL-8  IHC  images  from  E2-treated  and  Endo-R  MCF7L  xenograft  tumors.  (Scale  bars,  100  pm  and  50  pm,  respectively.)  (I)  Scatter  dot 
plots  of  IL-8  Allred  score  in  H.  Data  represent  means  ±  SEM,  *P  <  0.05,  **P  <  0.01,  ***P  <  0.001,  two-sided  t  test  for  indicated  comparisons.  (J)  Representative  IHC  images 
from  two  ER+  tumors  showing  low  (#1)  vs.  high  (#2)  FOXA1  and  the  negative  vs.  positive  IL-8  staining,  respectively.  (Scale  bar,  50  pm.)  (K)  Proportions  of  positive  vs. 
negative  IL-8  tumors  within  the  groups  of  tumors  showing  the  same  FOXA1  Allred  score  (AS).  Correlation  of  IL-8  positivity  and  FOXA1-AS  was  evaluated  by  Fisher's  exact  test. 
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pathways  (e.g.,  pAKT,  pMAPK,  and  pS6)  in  TamR  cells  was 
reduced  by  IL-8  knockdown  (Fig.  7C). 

To  further  investigate  the  relationship  of  IL-8  and  FOXA1  in 
endocrine  response,  we  established  a  series  of  inducible  MCF7L 
cell  lines  with  overexpression  of  YFP  (control)  or  FOXA1 
combined  with  concomitant  knockdown  of  luciferase  (control)  or 
IL-8  upon  induction  by  Dox.  As  a  result,  the  increased  IL-8  upon 
FOXA1  induction  was  substantially  reduced  by  coexpression  of 
IL8-shRNA  (Fig.  ID).  In  contrast,  FOXA1  induction  was  not 
altered  by  IL-8  knockdown  (Fig.  IE).  Without  Dox,  all  of  the 
MCF7L  stable  lines  showed  similar  sensitivity  to  endocrine 
treatment  (data  not  shown).  With  Dox,  IL-8  knockdown  alone 
increased  endocrine  sensitivity  in  the  ED  group  (Fig.  IF).  Con¬ 
versely,  FOXA1  overexpression  alone  decreased  the  endocrine 
sensitivity  to  all  of  the  antiestrogen  therapies.  Importantly,  the 
reduced  endocrine  sensitivity  by  overexpressing  FOXA1  could  be 
partially  reversed  by  concomitant  IL-8  knockdown,  suggesting  that 
IL-8  is  indeed  one  of  the  key  downstream  mediators  of  FOXA1  in 
conferring  endocrine  resistance.  Finally,  we  used  RPPA  to  mea¬ 
sure  the  signaling  changes  upon  concomitant  FOXA1  over¬ 
expression  and  IL-8  knockdown.  Among  the  proteins  up-regulated 
by  FOXA1  overexpression  in  MCF7L/FO}bkl  cells,  over  70% 
exhibited  reduced  expression  upon  concomitant  IL-8  knock¬ 
down  (Fig.  S8 B),  including  the  phosphorylated  proteins  of 
multiple  GFR  downstream  pathways,  such  as  AKT  (pAKT), 
JNK  (pc-Fos),  MAPK  (pMAPK),  JAK-STAT  (pSTAT3/6),  and 
mTOR  (pmTOR).  These  data  suggest  that  the  contribution  of 
IL-8  to  FOXAl-induced  endocrine  resistance  is  partially  me¬ 
diated  by  GFR  downstream  signaling  enhanced  by  high  FOXA1 
expression. 

Because  deregulated  IL-8  signaling  also  contributes  to  cancer 
cell  migration,  invasion,  and  metastasis  (40,  41),  we  next  evalu¬ 
ated  the  role  of  IL-8  in  cell  invasiveness.  We  found  that  IL-8 
knockdown  significantly  reduced  cell  invasion  in  MCF7L-TamR, 
but  not  MCFL-P,  cells,  which  are  much  less  invasive  at  baseline 
(Fig.  7G).  The  invasiveness  of  MCF7L-TamR  cells  was  also  re¬ 
duced  by  FOXA1  knockdown  (Fig.  S9A).  Both  600MPE-P  and 
TamR  cells  showed  stronger  invasiveness,  possibly  due  to  the 
constitutively  activated  RAS/RAF/MAPK  pathway.  IL-8  knock¬ 
down  partially  mitigated  the  invasiveness  of  both  600MPE-P  and 
TamR  cells  (Fig.  S9 B).  In  parallel,  FOXA1  overexpression  in 
both  MCF7L-P  and  600MPE-P  cells  enhanced  cell  invasion, 
which  was  abrogated  by  IL-8  knockdown  (Fig.  1H  and  Fig.  S9C). 
These  findings  support  a  role  for  IL-8  in  mediating  cell  invasion 
in  both  TamR  and  FOXAl-overexpressing  P  cells. 

Discussion 

In  characterizing  our  breast  cancer  Endo-R  cell  models  to  obtain 
clues  for  potential  mechanisms  of  endocrine  resistance  in  pa¬ 
tients,  we  discovered  gene  amplification  of  the  ER  pioneer  factor 
FOXA1  in  two  independently  derived  TamR  lines  of  MCF7 
cells,  and  we  found  FOXA1  overexpression  without  amplifica¬ 
tion  in  several  other  cell  lines  resistant  to  Tam  or  to  ED.  Recent 
studies  have  unveiled  gain-of-function  mutations  in  ESR1,  the 
gene  encoding  ER,  in  15-20%  of  metastatic  ER+  Endo-R  tu¬ 
mors  (42-45).  Genomic  amplification  or  overexpression  of 
FOXA1  may  be  another  mechanism  modulating  ER  activity  to 
promote  tumor  aggressiveness  and  endocrine  resistance.  We 
observed  FOXA1-CNG  and  amplification  in  20%  of  the  TCGA 
breast  tumors,  with  a  broader  FOXA1  -CNG  distribution  in  the 
luminal  B  subtype.  In  a  recent  study  reporting  genomic  profiling 
of  clinical  samples,  about  one-third  (7/20)  of  the  ER+  residual 
disease  showed  CN  changes  after  6  mo  of  neoadjuvant  anas- 
trozole  or  Ful  treatment  (46).  Interestingly,  compared  with  the 
baseline  tumors,  focal  amplicons  involving  the  FOXA1  or  ESR1 
gene  appeared  in  two  separate  cases  in  the  anastrozole  arm, 
supporting  clonal  selection  by  the  treatment  in  a  subgroup  of 
patients  as  a  mechanism  to  compensate  or  overcome  the  in¬ 
hibition  of  the  clinical  target/pathway.  These  data  provide  evi¬ 
dence  for  the  clinical  relevance  of  our  findings  in  the  Endo-R 
cell  line  models  and  further  suggest  that  the  genetic  alterations 
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Fig.  7.  IL-8  mediates  the  effect  of  FOXA1  on  cell  growth  and  invasion  in 
endocrine  resistance.  (A)  Cell  growth  within  5  d  in  MCF7L-P  and  TamR  cells 
with  IL-8  knockdown  by  two  different  sequences.  N.S.  knockdown  was  used 
as  normalization  control.  ( B )  A  stable  MCF7L-TamR/IL-8  cell  line  was  estab¬ 
lished  to  express  Dox-inducible  IL-8,  encoded  by  an  IL8  cDNA  without  3'-UTR 
sequence.  Two  different  IL8  siRNA  sequences,  targeting  either  the  IL8  coding 
DNA  sequence  (#1)  or  the  3'-UTR  region  (#2),  were  transiently  transfected 
into  MCF7L-TamR/IL-8  cells  ±  Dox  at  two  different  doses.  A  6-d  cell  growth 
measurement  was  performed  using  methylene  blue  staining.  Cell  growth 
under  N.S.  knockdown  was  used  as  the  normalization  control.  (C)  Western 
blots  of  GFR  downstream  signaling  mediators  in  MCF7L-TamR  cells  with 
siRNA  knockdown  of  N.S.  or  IL-8.  (D  and  E)  Measurement  of  IL-8  and  FOXA1 
mRNA  by  qRT-PCR  in  MCF7L  cell  lines  with  inducible  FOXA1  overexpression, 
or  concomitant  IL-8  knockdown  under  ±  Dox.  (F)  Cell  growth  within  7  d  in 
four  MCF7L  lines  with  Dox  induction  (0.5  ng/mL),  treated  with  E2  (as  control) 
or  antiestrogen  (ED,  Tam,  or  Ful).  Cell  growth  in  the  E2  group  was  set  as 
100%.  (G)  Cell  invasion  measurement  in  MCF7L-P  and  TamR  cells  transfected 
with  N.S.  or  IL-8-targeting  siRNAs.  Cells  were  seeded  onto  Matrigel-coated, 
24-well  Transwell  plates  and  cultured  for  48  h.  The  invading  cells  were 
counted  under  a  microscope  for  a  total  of  nine  random  fields.  Data  are 
presented  as  mean  number  of  cells  per  field.  ( H )  Cell  invasion  mea¬ 
surement  for  MCF7L/FOXA1  cells  ±  Dox  with  siRNA  knockdown  of  N.S. 
or  IL-8.  Cell  invasiveness  was  evaluated  as  above.  Data  represent  means  ± 
SEM,  *P  <  0.05,  **P  <  0.01,  ***p  <  0.001,  two-sided  t  test  for  indicated 
comparisons. 


in  the  ER  pathway  (e.g.,  FOXA1  and  ESR1)  might  drive  the 
outgrowth  of  rare  cell  populations  within  primary  tumors  that 
could  contribute  to  acquired  endocrine  resistance. 

In  addition  to  gene  amplification,  we  found  that  increased 
FOXA1  expression  occurred  at  the  mRNA  and  protein  levels  in 
other  Endo-R  cell  models  in  which  amplification  was  not  evi¬ 
dent.  The  epigenetic  and  posttranscriptional  regulation  of 
FOXA1  expression  in  breast  and  bladder  cancer  reported  by 
others  (47,  48)  might  also  apply  in  endocrine  resistance.  In  ad¬ 
dition,  a  recent  study  of  molecular  profiles  of  invasive  lobular 
carcinoma  identified  a  cluster  of  FOXA1  activating  mutations 
that  associated  with  its  expression  and  activity  in  promoting 
DNA  demethylation  of  its  binding  sites  (49).  The  chromatin 
binding  affinity  and  activity  of  FOXA1  can  also  be  modulated  by 
a  set  of  breast  cancer  risk-associated  SNPs  (50).  Collectively,  we 
speculate  that  there  are  multiple  mechanisms  by  which  FOXA1 
activity  can  be  up-regulated  in  the  setting  of  endocrine  resistance. 
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In  this  study,  we  showed  that  FOXA1  overexpression  in  ER+ 
breast  cancer  cells  activated  multiple  oncogenic  pathways, 
leading  to  endocrine  resistance  and  enhanced  cell  invasion. 
Conversely,  knockdown  of  FOXA1  in  TamR  cells  suppressed 
the  corresponding  oncogenic/GFR  downstream  signaling, 
leading  to  decreased  cell  growth  in  Endo-R  cell  lines.  Similarly, 
high  levels  of  FOXA1  have  been  shown  to  increase  the  growth 
of  prostate  cancer  cells  and  xenograft  tumors  and  to  correlate 
with  poor  prognosis  in  prostate  cancer  patients  (21,  51).  In 
breast  cancer,  high  levels  of  FOXA1  have  generally  been 
regarded  as  a  marker  of  good  prognosis  (15).  As  a  luminal  lineage 
determinant,  FOXA1  promotes  the  differentiation  of  normal 
mammary  epithelial  cells.  Likewise,  in  cancer  cells,  it  may  endorse 
a  classic  transcriptional  program  of  hormone  receptors  such  as 
ER,  resulting  in  a  more  differentiated  and  endocrine-sensitive 
phenotype. 

Previous  studies  from  our  group  and  others  have  shown  that 
during  ER+  disease  progression,  including  under  chronic  Tam 
treatment,  ER  switches  from  ligand  (E2)-dependent  to  ligand- 
independent  or  Tam-agonistic  signaling  and  a  transcriptional 
program  consistent  with  GFR  downstream  activation,  leading 
to  endocrine  resistance  (4,  5,  8,  37).  We  report  here  that  in¬ 
creased  levels  of  FOXA1  coordinate  at  least  partly  with  ER  in 
this  transcriptional  reprogramming,  leading  to  perturbed  gene 
signatures  and  signaling  pathways  associated  with  endocrine 
resistance.  As  such,  our  data  support  a  role  for  FOXA1  over¬ 
expression  in  more  aggressive  ER+  tumors,  which  is  in  line  with 
the  findings  of  high  levels  of  FOXA1  in  both  breast  and  pros¬ 
tate  cancer  metastases  (8,  52).  Moreover,  we  showed  that 
the  perturbed  genes  in  FOXA1 -overexpressing  MCF7L-P  cells 
were  enriched  for  the  predicted  genes  identified  by  FOXA1  and 
EGF-induced  ER  cistromes,  which  conforms  to  a  study  in 
MCF7  cells  showing  a  rapid  redistribution  of  ER  binding  me¬ 
diated  by  FOXA1  in  response  to  a  combination  of  mitogens  (8). 
As  such,  increased  levels  of  FOXA1  can  drive  ER  transcrip¬ 
tional  reprogramming  and  endocrine  resistance.  Strong  evi¬ 
dence  also  comes  from  prostate  cancer,  where  increased 
FOXA1  in  androgen-responsive  prostate  cancer  cells  facili¬ 
tates  AR-chromatin  binding  at  new  regions  and  promotes 
castration-resistant  and  androgen-independent  cell  growth 
(53).  Furthermore,  our  current  study  revealed  that  the  distinct 
cistrome  of  FOXA1  in  MCF7L-TamR  cells  was  enriched  for 
the  BCL11A  and  JUN/FOS  motifs.  Of  note,  we  have  pre¬ 
viously  shown  in  our  Endo-R  xenograft  tumors  the  increased 
activity  of  AP-1  (54),  the  transcription  factor  binding  to  the 
JUN/FOS  motif.  Further,  in  a  more  recent  study  integrating 
the  expression  data  from  our  Endo-R  xenograft  models  with 
the  previously  reported  growth  factor-dependent  ER  cistrome, 
and  using  functional  AP-1  blockade,  we  identified  AP-1  as  a 
key  determinant  of  endocrine  resistance  by  shifting  the  ER 
transcriptional  program  (55).  Altogether,  our  results  suggest 
complex  interplays  between  AP-1,  ER,  and  FOXA1  in  endo¬ 
crine  resistance  and  in  the  associated  genome-wide  transcrip¬ 
tional  reprogramming.  Importantly,  however,  in  our  Endo-R 
cell  models,  we  found  that  not  all  Endo-R  cells  with  increased 
FOXA1  maintained  ER  expression,  and  that  even  in  those 
maintaining  ER  not  all  of  the  transcriptional  reprogram  is 
influenced  by  knockdown  of  ER.  These  findings  suggest  that 
both  ER-dependent  and  ER-independent  (i.e.,  through  other 
transcription  factors)  mechanisms  may  underlie  the  impact  of 
increased  FOXA1  on  endocrine  resistance,  as  has  also  been 
suggested  in  prostate  cancer  androgen-deprivation  resistant 
models  (56).  In  the  context  of  this  study,  the  role  of  FOXA1  in 
mediating  AP- 1-dependent  gene  expression  in  an  ER-dependent 
or  -independent  manner  is  an  open  question  and  warrants 
further  study. 

How  exactly  does  FOXA1  at  high  expression  levels  induce 
endocrine  resistance?  Interestingly,  through  an  integrated  cis- 
tromic  and  transcriptomic  approach  and  functional  studies,  we 
identified  IL8  among  the  most  perturbed  genes  regulated  by 
FOXA1  in  an  ER-dependent  manner  in  TamR  cells.  Sub¬ 


stantial  evidence  indicates  that  increased  IL-8  levels,  through 
direct  effects  on  both  tumor  cells  and  tumor  microenviron¬ 
ment,  promote  survival  of  tumor-initiating  cells  (57),  tumor 
invasion  and  metastases  (40),  and  therapy  resistance  (41). 
However,  in  ER+  breast  cancer,  the  role  of  IL-8  remains  to  be 
determined.  It  has  been  reported  that  an  inflammatory  gene 
signature  identified  in  ER+  breast  tumors  is  associated  with 
poor  response  to  an  aromatase  inhibitor  (58).  We  find  that 
IL-8  mediates,  at  least  partially,  the  effect  of  increased  FOXA1 
on  cell  growth  and  invasion  in  our  Endo-R  cells.  IL-8  knock¬ 
down  effectively  inhibited  Endo-R  cell  growth  and  invasion, 
supporting  the  potential  of  IL-8  as  both  therapeutic  target  and 
biomarker  in  treating  Endo-R  tumors  with  high  levels  of  FOXA1 
and  IL-8. 

Collectively,  we  report  FOXA1  gene  amplification  and/or 
overexpression  in  Endo-R  cell  line  models.  Subclonal  evolution 
and  FOXA1/ER  transcriptional  reprogramming  may  coexist  as 
the  underlying  mechanism  of  endocrine  resistance.  IL-8  signaling 
is  one  of  the  components  embedded  in  the  FOXA1/ER  tran¬ 
scriptional  reprogramming  and  provides  a  potential  therapeutic 
target  for  ER+  tumors  with  increased  FOXA1. 

Materials  and  Methods 

The  Endo-R  derivatives  were  developed  from  P  cells  of  MCF7L  from 
M.  Lippman  (Sylvester  Comprehensive  Cancer  Center,  Miami,  FL),  >- 

600MPE  (J.W.G.),  ZR75-1  (American  Type  Culture  Collection),  and  BT474  o 

(AstraZeneca),  using  the  method  we  previously  reported  (59).  The  MCF7RN  2 

Endo-R  cell  model  was  kindly  provided  by  R.  Nicholson  and  J.  Gee,  Cardiff  j 

University,  Cardiff,  UK.  All  of  the  cells  were  authenticated  and  the  P  cells  S 

were  maintained  in  RPMI/1640  (MCF7,  ZR75-1)  or  DMEM/high-glucose 
(600MPE,  BT474),  supplemented  with  10%  (vol/vol)  heat-inactivated  FBS 
and  1%  (vol/vol)  penicillin/streptomycin/glutamine  (PSG).  The  Endo-R 
cells  were  kept  in  phenol-red  free  (PRF)  medium  supplemented  with  10% 
(vol/vol)  heat-inactivated  charcoal-stripped  (CS)-FBS  and  1%  (vol/vol) 

PSG,  with  (for  TamR)  or  without  (for  EDR)  the  addition  of  100  nM  4-OH- 
Tam  (H7904;  Sigma).  The  Dox-inducible  FOXA1 -overexpressing  cell  lines 
were  established  using  a  lentiviral  cDNA  delivery  system  from  X.  Pan, 
(Novartis,  Cambridge,  MA)  and  maintained  by  200  |ag/mL  Geneticin 
(Invitrogen).  The  Dox-inducible  shlL-8  knockdown  cell  lines  were  estab¬ 
lished  using  the  pINDUCER  system  (60).  Exome-seq  and  CN  analysis,  FISH 
and  gPCR  assay.  Western  blotting,  animal  studies,  RNA  interference,  cell 
growth  assay,  RNA-seq  and  gene  expression  analysis,  Kaplan-Meier 
curves,  RPPA  and  signaling  pathway  analysis,  integrated  ChIP-seq  and 
RNA-seq  data  analysis,  qRT-PCR,  integrative  cistromes  analysis,  ChlP- 
qPCR,  ELISA,  IHC,  immunofluorescence  staining,  and  cell  invasion  assay 
are  described  in  SI  Materials  and  Methods.  Statistical  analysis  of  in  vitro 
assays  was  based  on  at  least  triplicated  data  using  R  software  (v2.1 3.0)  or 
GraphPad  Prism  (v5.04).  All  experiments  were  repeated  at  least  three 
times.  Quantitative  data  are  shown  as  mean  ±  SEM  from  triplicates  or 
quadruplicates.  Significant  difference  (P  <  0.05)  was  determined  by 
ANOVA  or  Bonferroni  post  hoc  tests  (multiple  testing  corrected). 

Animal  care  and  animal  experiments  from  this  study  were  in  accordance 
with  and  approved  by  the  Baylor  College  of  Medicine  Institutional  Animal 
Care  and  Use  Committee. 
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Background:  Metastasis  in  ER-positive  (+)  breast  cancer  (BC)  occurring  years  to  decades  after  initial  diagnosis 
presents  a  daunting  challenge  for  clinical  care  and  preclinical  research  due  to  limited  known  key  players  and 
experimental  models.  FOXA1  is  a  pioneer  factor  for  ER-chromatin  binding  and  function,  and  is  highly 
expressed  in  ER+  BC  metastases,  yet  the  underlying  mechanism  is  unclear.  Tumor-secreted  proteins  play  a 
crucial  role  in  the  reciprocal  interplay  between  cancer  cells  and  host  microenvironmental  factors  at  both  primary 
and  secondary  sites.  We  hypothesized  that  high  FOXA1  provokes  an  ER-dependent  transcriptional  program  that 
includes  a  unique  pro-tumorigenic  secretome  essential  for  promoting  ER+  BC  metastasis.  Methods:  A  lentiviral 
doxycycline  (Dox) -inducible  FOXA1  overexpression  vector  and  a  dual  luciferase/GFP  (LG)  tracking  vector 
were  integrated  to  construct  a  stable  MCF7-LG/FOXA1  cell  model.  Ovariectomized  nude  mice  bearing  MCF7- 
LG/FOXA1  xenografts  in  the  presence  of  exogenous  estrogen  (E2)  were  randomized  to  ±  Dox,  each  with 
continued  E2,  E2  deprivation  (ED),  or  tamoxifen  (Tam).  Survival  surgery  removing  the  therapy-naive  (E2  arm) 
and  relapsed  (ED/Tam  arms)  tumors  was  performed  when  tumors  reached  -1000  mm3.  All  mice  then  received 
ED/Tam  ‘adjuvant’  therapy,  with  longitudinal  luminescence  imaging  to  monitor  local/distant  recurrences.  Mice 
were  or  will  be  euthanized  at  the  ethical  end-point.  Integrative  bioinformatics  was  performed  using  RNA-seq 
and  FOXA1/ER  ChIP-seq  data  from  our  preclinical  models  to  identify  secretome  targets  for  functional 
intervention.  Times  to  tumor  regression  (TTR)  and  progression  (TTP)  were  defined  by  when  the  tumor  reached 
half  or  twice  the  volume  at  randomization.  Results:  Median  (m)  TTR  was  achieved  in  ED  (31/34  days,  -/+Dox, 
P  =  0.184)  but  not  in  Tam  groups  —  Tam  delayed  tumor  growth  but  failed  to  prevent  progression  in  all  mice 
with  mTTP  of  94/93  days  (-/+Dox,  P  =  0.517).  Despite  no  difference  in  mTTP  at  Tam-/+Dox,  a  quarter  of 
+Dox  tumors  (3/12)  had  volume  doubled  by  day  11.  No  metastases  were  observed  by  imaging  in  any  of  the 
mice  before  surgery  (‘neoadjuvant’  setting).  Local  relapse  and  lymph-node/lung  metastases  were  detected  after 
surgery  (‘adjuvant’  setting).  At  day  90  in  the  adjuvant  Tam  group  with  previously  relapsed  tumors,  +Dox  mice 
succumbed  to  metastasis  more  often  than  -Dox  mice  (7/8  vs.  3/10,  P  =  0.023).  Compared  to  the  adjuvant 
Tam+Dox  mice  with  previous  therapy-naive  tumors,  the  Tam+Dox  with  previously  relapsed  tumors  showed 
higher  distant  metastasis  rate  (7/8  vs.  5/14,  P  =  0.026).  Analysis  of  the  ED  setting  is  pending  due  to  late 
recurrence.  Data  integration  and  functional  study  revealed  a  set  of  cytokines,  growth  factors,  and  extracellular 
matrix  components  (including  IL-8,  CTGF,  and  LOX),  regulated  by  FOXA1  often  in  conjunction  with  ER,  that 
are  highly  involved  in  FOXAl-induced  metastasis.  Global  secretome  profiling  by  mass  spectrometry  and  target 
validation  are  ongoing.  Conclusions:  FOXA1  overexpression  increases  metastatic  potential  in  ER+  BC.  We 
established  a  pertinent  metastatic  xenograft  mouse  model  to  characterize  a  pro-metastatic  secretome  with 
diagnostic  and  therapeutic  potential  for  treating  metastatic  ER+  BC. 


